Power Law in Hadron Production 



Marek GazdzickiQ 



CERN, Geneva, Switzerland 
and 

Institut fur Kernphysik, Universitat Frankfurt 
D-60486 Frankfurt, Germany 



Mark I. Gorensteirj^ 

Institut fur Theoretische Physik, Universitat Frankfurt 
D-60054 Frankfurt, Germany 
and 

Bogolyubov Institute for Theoretical Physics, Kiev, Ukraine 



In high energy p(p)+p interactions the mean multiplicity and transverse mass spectra 
of neutral mesons from rj to T (m = 0.5 -f- 10 GeV/c 2 ) and the transverse mass spectra 
of pions (m,T > 1 GeV/c 2 ) reveal a remarkable behaviour: they follow, over more than 
10 orders of magnitude, the power-law function: C ■ m^. The parameters C and P 
are energy dependent, but similar for all mesons produced at the same collision energy. 
This scaling resembles that expected in the statistical description of hadron production: 
the parameter P plays the role of a temperature and the normalisation constant C is 
analogous to the system volume. The fundamental difference is, however, in the form of 
the distribution function. In order to reproduce the experimental results and preserve 
the basic structure of the statistical approach the Boltzmann factor e~ E ' T appearing in 
standard statistical mechanics has to be substituted by a power-law factor (E* /A)~ p . 
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It is well established that basic properties of hadron production in high energy collisions approximately 
follow simple rules of statistical mechanics. A principal assumption of statistical models of strong inter- 
actions jjj states that the single particle energy distribution in the local rest frame of hadronizing matter 
follows a Boltzmann form: 

dn ( E*\ . . 

exp -— , (B.l) 



dp* \ T 



where E* = \J to 2 + p* 2 = my • cosh(y*) is the hadron energy (% = \J m 2 + pj, 
and y* is the local hadron rapidity) and T is a common (for all particles) temperature parameter usually 
extracted by a comparison with the experimental data. In the fixed rest frame (e.g. in the cm. frame 
of the colliding particles) the energy E* equals to m T ■ cosh(y — y'), where y and y' are the rapidities of 
considered particle and the 'hadron fluid elements' (fireballs), respectively. 

Integration of hadron distribution along the collision axis (i.e. over y and y' with an arbitrary rapidity 
distribution function, f(y'), of fireballs) leads to an approximately (for uit » T) exponential form of 
transverse mass spectrum jlj : 

exp - . (B.2) 
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According to this approach the final state hadrons should obey this distribution providing that collective 
transverse motion of hadronizing matter and contributions from resonance decays are neglected. Further 
integration over transverse mass yields an expression for the mean hadron multiplicity, N(m), which is 
also (for to >> T) governed by the Boltzmann factor: 

(in \ 



The exponential distributions (B.2, B.3| ) are confirmed in high energy particle collisions by numerous 



experimental results on px spectra in the transverse momentum region pt < 2 GeV/c and on hadron 
yields at m < 2 GeV/c 2 j!| (a correction for contributions from resonance decays is necessary for a det ailed 
comparison with data). For a given hadron species, h, the normalisation factors in Eqs. (]B.2| , |B.3| ) are 
proportional to a volume parameter V (a sum of the proper volume of 'hadron fluid elements'), which 
has to be common for all hadrons, following the assumption of statistical production. Additional factors 
are degeneracy factor g = (2j + 1), where j is the particle spin, and a chemical factor exp(/x/j/T), 
which accounts for material conservation laws in grand canonical approximation. The chemical potential 
Hh equals to zero for neutral hadrons (i.e. for hadrons with zero valu es of the conserved charges). 



Consequently, for neutral hadrons the exponential distributions (B.2, |B.3| ) with a common temperature 
parameter and a common volume parameter (multiplied by the degeneracy factors g) are expected to be 
valid. 

However, in high energy p+p interactions the measurements of ir° meson spectra in a high pt {pr > 
2 GeV/c) domain show a strong deviation from the exponential shape. In this region the distribution 
follows a power-law dependence J 31] . In the standard, QCD-based, approach this behaviour is attributed 
to the hard parton scattering Mj. A strong violation of the Boltzmann behaviour in the high mass 
region is also seen in the case of mean hadron multiplicities. The hadron yields from the ?7-meson 
(to = 0.55 GeV/c 2 ) to the J/-0-meson (to = 3.1 GeV/c 2 ) follow @] approximately an exponential law 



(B.3). This is, however, not true any more for the T meson (to = 9.5 GeV/c ). The experimental ratio 



of T to J/ip mean multiplicity at 800 GeV/c is |( 

(T) 



whereas Boltzmann's exponential law predicts: 



10" 
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for typical value of T = 0.170 GeV. Thus the standard statistical model underpredicts the experimental 
data by more than 10 orders of magnitude! 



In this letter we discuss a possibility to extend the statistical approach of hadron production in el- 
ementary hadronic interactions to the high (transverse) mass domain. Based on experimental data we 
postulate that the energy spectrum in the local rest frame of hadronizing matter is given by a power-law 
distribution: 



dn (E* 
dp* ~ V A 



(P.l) 



instead of the exponential Boltzmann distribution (B.l). We assume further that the remaining structure 



of the statistical approach is unchanged. Two new parameters appear in the proposed statistical power- 
law model: a scale parameter A and an exponent P, both are assumed to be common for all hadrons. 
We do not attempt here to introduce proper canonical treatment of material conservation laws needed 
for description of charged hadrons in small systems J?],||]. Therefore, we limit our analysis to the case of 
neutral mesons only. Further on, in order to be able to neglect effects of large resonance widths in the 
statistical treatment we consider only narrow mesonic states. 



Integration of Eq. (P.l) over longitudinal motion results in a power-law transverse mass distribution: 

p(m T ) = ^ = C ■ m T P , (P.2) 

where C is a normalisation parameter in which we absorbed the dependence on the scale (A) and volume 
(V) parameters. As both A and V are assumed to be common for all hadrons, so does the normalisation 
parameter C . 

The integration of Eq. ( |P.2| ) over transverse mass yields: 

p(m) = — ^| ■ N(m) = C ■ m- p , (P.3) 
g to 

where N(m) is the hadron multiplicity. 



Several nontrivial predictions of the statistical power-law model follow from Eqs. (P.2) and (P.3) 



• both transverse mass spectra and hadron yields should obey a power-law behaviour; the power P 



in the mT _ distribution, p(nvr) (P-2), and in the hadron yield formula, p(m) (P.3), should be the 
same, 



• the normalisati on co nstant C also should be the same for both TOy-distribution (P.2) and hadron 



yield spectrum (P.3) 



• the power P and the normalisation constant C should be universal (equivalent to the temperature 
and volume parameters in the standard statistical approach), i.e. they are the same for different 
hadron species. 

In order to test these predictions we plot in Fig. 1 experimental results on m^-spectra and hadron 
yields for p+p interactions at y/s — 27 4- 30.8 GeV as a function of tot and m, respectively. Within 
the statistical power-law model, results on both p(m) vs m and p(nvr) vs tot should follow the same 
dependence: C ■ rn7 T -, . 

Full dots in Fig. 1 indicate results on p(m) calculated for the mean multiplicity of hadrons: 
t],uj,4> pi, J/if), ip' and T ||, whereas triangles indicate the data on TOr-spectra, p(mr), of neutral 
pions fiOff T The measurements of J/ip,ip' and T || were performed for p+S interactions at 800 GeV/c. 
The extrapolation to p+N interactions was done assuming the A 092 dependence on nuclear mass num- 
ber. The measured energy dependence of midrapidity yield || was used to calculate multiplicities at 
•\/s w 30 GeV. The ir° measurements are done at central rapidities only. In order to calculate the rapid- 
ity integrated distributions needed for comparison with the model, we use the following approximation: 
dN/drriT — d 2 N / (drriTdy) ■ Ay, where Ay = 1. Measurements of the transverse momentum spectra 
of 7T° mesons performed by several experiments at the ISR [[h| differ typically by a factor of about 2. 
This together with similar uncertainty introduced by the extrapolation procedure for rapidity integrated 
distributions, yield a relatively large systematic error on the absolute normalisation of ir° spectra. An 
additional bias of a similar magnitude is possible due to a missing correction for contributions from res- 
onance decays. A more detailed analysis is obviously needed in the future. Nevertheless, we observe a 



surprising agreement over more than 10 orders of magnitude of the experimental results with the expecta- 
tions following from the statistical power-law model. Typical deviations of the experimental points from 
the universal power-law dependence are of about factor 2, which is a similar magnitude as the systematic 
errors estimated above. In the power-law fits to the data the statistical and 50% systematic errors were 
added in squares. The values of the parameter P (= 10) and the noramlization constant C resulting from 
the separate fits to the mj- spectra and hadron multiplicities are similar (for details see Table 1). The 
solid line in Fig. 1 indicates a fit to the mult iplicity d ata. 

As seen from Fig. 1 the power-law model (P.2, P.2) works rea sonably well even for intermediate values 
of ni(T) (0.5 -j- 2 Ge V/c 2 ), where the exponential distributions (B.2, B.3|) are normally assumed, e.g. the 
power-law function ( [P.3] ) seems to describe also the data for r\ and (f> multiplicities. This fact requires a 
further study. 

In order to check whether the observed power-law ray-scaling is valid also for higher collision energies 
we compiled data on J/ip, ip' jyj, T(ls), T(2s) and T(3s) Jl2| in p+p interactions at y/s — 1800 GeV 
(the highest energy accessible today) . The pr distributions for the quarkonia are measured at midrapidity 
and they allow to calculate p(m*r) spectra in this region. The resulting distributions are indicated by 
open symbols in Fig. 1 together with the line showing the power-law fit (P = 8, for details see Table 1) 
to these data. It is seen that the distributions of quarkonia produced in p+p interactions at y^s — 1800 
GeV follow approximately the same power-law function. The relative comparison of data for different 
quarkonia may be biased by a systematic error of about 50% due to the use of midrapidity spectra instead 
of spectra integrated over all rapidities. Therefore, as previously, in the fit the square of the error was 
calculated as a sum of the squares of statistical and systematic (50%) errors. 

There are no data on 7r° meson production in p+p interactions at y/s = 1800 GeV. However, there are 
measurements of ir° spectra performed close to midrapidity (y* w 1.4) at ^/s = 540 GeV which extend 
up to pr = 40 GeV/c The resulting p(rriT) distribution is shown in Fig. 2. The fit of the power-law 
function indicated by a solid line yields P = 8 (for details see Table 1), the value similar to one obtained 
for quarkonium spectra at 1800 GeV. The power-law fits to 30 and 1800 GeV data, indicated in Fig. 
2 for a comparison, are both below the spectrum at 540 GeV. This observation may suggest a possible 
violation of the scaling for ir° mesons at 1800 GeV. This is because one may expect a monotonic increase 
of p(txit) at fixed m-r with increasing collisions energy and this expectation is not obeyed by the 540 GeV 
data. We note, however, that the measurements of it spectrum at high pt region were indirect due to 
difficulty to resolve single photons from ir° decay. Therefore, a direct measurements of tt° spectrum at 
1800 GeV are necessary in order to verify a validity of the scaling in this domain. 



What is the origin of the power-law ray-scaling observed in p(p)+p interactions and not (yet) seen in 
the collisions of heavy nuclei JlJ]? The production of heavy hadrons (like J/ip and T mesons) as well as 
the power-law behaviour of pion spectra at high px are fitted by QCD inspired models [EUTq] . However 
different assumptions and parameters enter in these models for ir° and quarkonia description. Thus the 
observed power-law ray-scaling seems to be an unexpected and interesting feature of the data which 
requires explanation. In this context we list known to us developments which can be helpful in attempts 
to understand the power-law ray-scaling. It was recently pointed out Jl(f| that quantum fluctuations of 
the string tension may account for the thermal features of the hadron spectra. A possible appearance of 
the power-law instead of Boltzmann spectrum was suggested within thermal field theory |17| in the low 
temperature and high mass limit. Generalization of standard thermodynamics to non-extensive systems 
|l8"f leads to distribution of energy in the canonical ensemble with the power-law tail for high energies. 

Finally we repeat the basic results presented in this letter. Properly normalized multiplicities and ray 
spectra of neutral mesons from 77 to T (m = 0.5 -7- 10 GeV/c 2 ) and my spectra of ir° meson (mj- > 
1 GeV/c 2 ) obey the ray-scaling. The scaling function has approximately the power-law form: m^. 
The values of parameter P and normalisation constant C resulting from the fits to data on production 
of different mesons at fixed collision energy are similar. This scaling behaviour resembles that expected 
in statistical mechanics: the parameter P plays the role of temperature and the normalisation constant 
C is analogous to the system volume. Thus the basic modification of the statistical approach needed to 
reproduce the experimental results on hadron production in p(p)+p interactions in the large TO(y) domain 
is the change of the shape of the distribution function. The Boltzmann function e~ E l T appearing in the 
standard statistical mechanics has to be substituted by a power-law function (E* /A)~ p . 
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Table 1 The values of the parameters (P and C) resulting from the power-law fits (p(m^) = C-m,^) 
to various data sets discussed in this letter (for details see text). The values of x 2 /ndf are also given. 



data set 


P 


C-10 2 (GeV<P-3)) 


X 2 /ndf 


neutral meson multiplicity at =30 GeV 


10.1 ± 0.3 


3.6 ± 1.0 


1.02 


mT-spectra of n° at =30 GeV 


9.8 ± 0.1 


6.1 ± 0.9 


1.4 


m-T-spectra of quarkonia at 1800 GeV 


7.7 ± 0.4 


3.4 (+4.4,-2.0) 


1.13 


m-T-spectra of tt° at 540 GeV 


8.1 ± 0.1 


129 ± 15 


1.3 
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FIG. 1. Mean multiplicity of neutral me sons (full dots) scaled to p(m) (P.3) and transverse mass spectra of 



7r" mesons (full triangles) scaled to pimr) ( |P.2| ) produced in p+p interactions at y^i w 30 GeV as well as the 
quarkonium p(rriT) spectra (J/ip - crosses, ip - diamonds, T(ls) - stars, T(2s) - triangles, T(3s) - squares) for 
p+p at t/s — 1800 GeV. The solid lines indicate power-law fits to the data. 
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FIG. 2. The p(mr) spectra of 7r° meson for p+p interactions at yfs = 540 GeV. The solid line indicates a 
power-law fit to the 540 GeV data, whereas the two dashed lines show the corresponding fits to the results at 30 
and 1800 GeV (see Fig. 1). 
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